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Abstract 

We use the shell model Monte Carlo method to calculate complete 0/lp- 
shell response functions for Gamow- Teller (GT) operators and obtain the 
corresponding strength distributions using a Maximum Entropy technique. 
The approach is validated against direct diagonalization for ^^Ti. Calculated 
GT strength distributions agree well with data from and (p, n) reactions 

for nuclei with ^ = 48 — 64. We also calculate the temperature evolution 
of the GT+ distributions for representative nuclei and find that the GT-(_ 
distributions broaden and the centroids shift to lower energies with increasing 
temperature. 

PACS numbers: 21.60.Cs, 21.60.Ka, 27.40.+Z, 23.40.-s 
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I. INTRODUCTION 



The Gamow- Teller (GT) properties of nuclei in the medium mass region of the periodic 
table are crucial determinants of the precollapse evolution of a supernova |jl|] . The core of a 
massive star at the end of hydrostatic burning is stabilized by electron degeneracy pressure 
as long as its mass does not exceed the appropriate Chandrasekhar mass Mqh- If the core 
mass exceeds Mch, electrons are captured by nuclei. For many of the nuclei that determine 
the electron capture rate in this early stage of the presupernova Gamow- Teller (GT) 
transitions contribute significantly. Due to insufficient experimental information, the GT_|_ 
transition rates have so far been treated only qualitatively in collapse simulations, assuming 
the GT_|_ strength to reside in a single resonance whose energy relative to the daughter ground 
state has been parametrized phenomenologically 0; the total GT_|_ strength has been taken 
from the single-particle model. However, recent {n,p) experiments show that the 

GT_(_ strength is fragmented over many states, and that the total strength is significantly 
quenched compared to the single-particle model. (A recent update of the GT+ rates for use 
in supernova simulations assumed a constant quenching factor of 2 0.) 

In this paper, we describe our calculations of Gamow- Teller strength distributions in iron 
region nuclei: the shell model Monte Carlo (SMMC) technique is used to obtain the response 
functions of the Gamow- Teller operators in the full Ohu /p-shell model space. These response 
functions are related to the strength distributions through an inverse Laplace transformation, 
which we carry out using a Maximum Entropy method. 

Our starting point is the interacting shell model |^, which gives an accurate and con- 



sistent description of the properties of light nuclei |]10|,|ri| when an appropriate interaction 
is used. In the shell model, nucleons occupy a spectrum of single-particle orbitals that are 
formed by the presence of an assumed mean field. These nucleons interact through a resid- 
ual effective interaction, which is derived from a realistic nucleon-nucleon potential through 
the G-matrix formalism [|12|. The resultant interaction matrix elements require some min- 
imal tuning to optimally account for known spectroscopic properties. In the conventional 
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approach, the solution to the shell model is obtained by diagonalizing the nuclear Hamil- 
tonian in a suitably chosen basis of many-particle configurations. Since the Hamiltonian 
matrix to be diagonalized grows combinatorially with the size of the single-particle basis 
and the number of valence nucleons, realistic calculations are feasible in the full /p-shell 
only for nuclei with A < 50. Hence, the traditional calculation of various nuclear properties 
for medium-heavy and heavy nuclei lies beyond the scope of direct-diagonalization methods 
except in a severely truncated model space. 

The SMMC method [[TsHlBI scales more gently with the problem size than do traditional 
direct-diagonalization techniques, allowing larger, and hence more realistic, calculations. 
This method exploits the fact that most of the billions of configurations in nuclei are unim- 
portant for general nuclear properties, so that only a subset of the relevant configurations 
needs to be sampled. Observables are calculated as thermal averages in a canonical ensemble 
of nuclear configurations, so that nuclei at finite temperature can be studied quite naturally. 

SMMC methods were used in the first complete Ohu calculations for a number of ground- 
state |T7|-P^, and finite-temperature properties of mid- fp shell nuclei. These studies 
used both the Richter-Brown [^] and the KB3 residual interactions. For the purposes of 



investigating Gamow- Teller transitions, the KB3 interaction (obtained by minimally modi- 
fying the monopole strength in the original Kuo-Brown matrix elements [^) is well-suited 
for full Ohuj studies throughout the lower- /p shell region |j23|. Observables that have been 
calculated with this interaction in the SMMC approach include the energy (H), the total 
B{E2), B{M1), GT strengths, and various pairing properties; the calculated ground-state 
properties compare very well with experiment. Importantly, these studies showed that the 
experimentally observed quenching of the total GT strength is consistently reproduced by 
the correlations within the full /p-shell if a renormalization of the spin operator by the fac- 
tor 0.8 is invoked [^,^. The same renormalization factor had already been deduced from 



sd-shell [11 1 and /p-shell nuclei with A < 49 [25,26| and thus appears to be universal 



In Section II, we reveiw the SMMC method and its application to response functions. 
We apply a Maximum Entropy (MaxEnt) method to perform the required inverse Laplace 



transform of the SMMC response functions; our implementation of MaxEnt for SMMC is 
discussed in Section III. Section IV includes a validation of these methods against direct 
diagonalization for GT transitions in ^^Ti, and we present GT strength functions for several 
heavier nuclei in the /p-shell {A = 48 — 64) where experimental data are available. We also 
discuss the evolution of these distributions with temperature. A brief conclusion follows in 
Section V. 



II. THE SHELL MODEL MONTE CARLO METHOD 

The SMMC method is based on a statistical formulation of the nuclear many-body 
problem. In the finite-temperature version of this approach, an observable is calculated as 
the canonical expectation value of a corresponding operator ^ at a given temperature T 
and is given by [p!^-p!6| 



where U = exp(— is the imaginary-time many-body propagator, TtaU is the canonical 
partition function for A nucleons, H is the shell model Hamiltonian, and (3 = 1/T is the 
inverse temperature. 

In terms of a spectral expansion, the total strength of a transition operator A is then 
given by the following expectation value: 



where \i) (|/)) are the many-body states of the initial (final) nucleus with energy Ei (Ef). 
The total strength from the ground state can be obtained by choosing a sufficiently large 
value for (3 such that only the ground state contributes due to the Boltzmann weight. 

In addition to the "static" strength [Eq. (2)], one can calculate for an imaginary-time r, 
the response function, Ra{'t)i which describes dynamical behavior and contains information 
about the nuclear spectrum: 
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RAr) = [AHr)Am) = = ^-^^ . (3) 



The strength distribution 



c E^fS{E-Ef + E,)e-^^^\{f\A\^)\' 

^A[E) = (4) 



is related to -R^(r) by a Laplace Transform, 

/oo 
SAE)e--''dE. (5) 
-oo 

Note from Eq. (3) that ground-state to ground-state transitions require large (/5 — r) 
in addition to large (3. The large-r behavior of i?^ allows, in principle, a measurement of 
the specific transition between the ground state and the lowest allowed final state by the 
operator; the slope of loge[-R(T)] in this limit provides the transition energy, and the intercept 
measures the transition strength. 

The SMMC canonical expectation values are based on the discretization of the many- 
body propagator, e~^^ , into a finite number of "time" slices, Nt, each of duration 
A/5 = (3/Nt. At each time slice the many-body propagator is linearized via the Hubbard- 
Stratonovich transformation p7|j28[] ; observables are thus expressed as path integrals of 
one-body propagators in fluctuating auxiliary fields. The integration is carried out by a 
Metropolis random walk p9| . 

To circumvent the "sign problem" encountered in the SMMC calculations with realistic 



interactions, we use the extrapolation procedure outlined in Refs. Yet another, but 

distinct, source of the sign problem is an odd number of nucleons in the canonical expectation 
values |jl6|. We overcome this problem by a number-projection technique, first employed in 



T8| and subsequently used in ||T6[, that allows us to extract information concerning odd-A 



nuclei from the neighboring even-even system. 



III. THE METHOD OF MAXIMUM ENTROPY 



Once we have the Gamow- Teller response functions, they must be inverted to obtain 
strength distributions. The inverse of the Laplace transform (5) required to extract the 
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strength functions is an ill-conditioned numerical problem The kernel (which in this 



case is e ) acts as a smoothing operator and thus the solution, for which the kernel must 
be inverted, will be extremely sensitive to small changes, {i.e., to errors) in the input data. 



In this section, we describe a Maximum Entropy procedure to carry out the inversion |T6 

Consider the x^-deviation of the data, = i?(r = iA(3), with errors, a,, from the fit 
values Fi{S} produced by the trial inverse and obtained according to Eq. (H), 



x'{s} = E 



,n-F,{S}. 



(6) 



Direct minimization of is numerically stable only in the simplest of circumstances (such 
as few-parameter data fitting). Combining with some other auxiliary well-conditioned 
functional, P {S}, such that P{S} has a minimum at the smooth solution, S{E), and 
penalizes strongly oscillating functions, leads to a compromise between fitting the data and 
the expected smoothness of the inverse. Thus one minimizes the joint functional 

1 



-X'{S} + P{S} 



(7) 



The functional P{S} is chosen as the information theoretic entropy, 
P{S} = aJ dE 



m{E)-S{E) + SiE) In (^^) 



(8) 



where m{E) is a default model and a is an adjustable parameter that both specify the a 
priori knowledge of S{E). 

In order to minimize the functional (0), we employ the technique of Ref. which 
involves an iterative sequence of linear programming problems. We first expand Eq. (|^) to 
second order in S{E) about some positive function f{E) to obtain 



P{f\S} = aJdEUm-^) + 



In 



m 



1 



S 



2/ 



(9) 



If the true minimum S{E) of the non-quadratic functional in Eq. (^) is taken as a point of 
expansion of f{E) in [Eq. (^], then it also gives the minimum of the corresponding quadratic 
functional 
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S{E) = min 



^x'{a} + P{S\a} 



(10) 



Since we require extraction of positive strength function, we iterate while retaining par- 
tially the result of the previous iteration as, 



5(^^+1) = min 

5>0 



(11) 



'lx'{S} + p{f''Us}\ , 

with 

/(")(^) = ^S^^-'\E) + (1 - OS^'^KE) , (12) 

and the default model as the starting approximation to S, 

S^^'^E) ^ S^-^\E) = m{E) . (13) 

The rate of convergence and stability are controlled by the mixing parameter < ^ < 1; 
a value of ^ = 0.3 is a reasonable choice to guarantee stabihty. Typically, convergence to 
the "true" solution is obtained in less than 40 iterations. In this way, the minimization 
of a general functional that is intrinsic to a Maximum Entropy approach is reduced to an 
iterative procedure in which each step requires the minimization of a quadratic functional 
with linear inequality constraints. 

Some general remarks regarding this inversion technique are called for. Since R{t) is 
calculated at discrete values of imaginary time and, in principle, up to an imaginary time 
/3, the smallest energy that can be resolved in S{E) is of order 1//3, and the largest is the 
inverse of the discretization size, l/A/3. In practice, numerical noise forces a cut-off in the 
largest r value that can be used, thus decreasing the energy resolution. 

As we mentioned above, the default model can be chosen by investigating the charac- 
teristics of the response function. From Eq. (3), one sees that dln[R{T)]/ dT\r=o gives the 
centroid of the distribution in the parent nucleus, and thus in the case of the GT+ operator 
we choose for the default model a Gaussian with a peak at this energy and with a width of 
1.5 — 2 MeV; this width can be estimated from (f ln[R{T)]/ dT^\T=o- The parameter a is the 



inverse of the total strength of the distribution, and is calculated from the default model as 
a = [J dEm{E)]'^ . In the case of the GT_ operator, we make a better guess for the default 
model by including some features of the distribution. Experimental distributions typically 
have three regions: the T = T^, and T = + 1 regions distributed around 6 MeV and 
12 MeV, respectively, and a more fragmented region at lower energies. We choose for our 
GT_ default model two gaussians with the same widths, each centered at the appropriate 
energy. The lower energy part of the distributions is governed by the high r region of the 
response function. Although this region of the response function is sometimes contaminated 
by large statistical fluctuations, the reconstruction tends to give a low-energy peak that well 
describes these more discrete transitions. 



IV. GAMOW-TELLER STRENGTH DISTRIBUTIONS 

The GT operators are defined as GT± = J2i <^i't^^ where <Ji is the Pauli spin operator for 
nucleon / and {ti') is the isospin lowering (raising) operator that changes a neutron (pro- 
ton) into a proton (neutron); they thus describe charge- changing decay modes. GT strength 
distributions play an important role in two very different contexts. In the astrophysical 
context, medium- heavy nuclei at a finite temperature in the core of a pre-supernova capture 
electrons. A strong phase space dependence makes the relevant electron capture rates more 



sensitive to GT distributions than to total strengths [0,0 and thus necessitates complete 
{]huj calculations of these distributions. GT strengths are also important in studies of double 
beta decay @]. The two-neutrino mode of this decay, which provides important confidence 



in extracting the neutrino mass from zero-neutrino decay experiments, is equivalent to a 
description of the GT strength functions from the ground states of the parent and daughter 
nuclei. Thus, any reliable calculation of the two-neutrino matrix element must accurately 
describe these strength distributions. 

In the following sections we demonstrate and validate the MaxEnt method for the GT 
operator by comparing our results with direct diagonalization. We then compare our results 
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with experimentally obtained distributions for various /j9-sliell nuclei. In what follows we 
will use the renormalized GT operator corresponding to GT-I-/1.26 ||I9| , |23 |. 



A. Comparison with direct diagonalization 

Direct-diagonalization results in the complete /p-shell can be obtained for nuclei with 
A < 48. We choose ^^Ti for a comparison and in Fig. 1, we show our results for this nucleus. 
The lower left-most panel shows the GT+ response function, R{t), for ^^Ti as measured in 
the parent and the middle lower panel shows the extracted strength distribution, S{E), in the 
daughter, ^^Sc. Also shown in the same panel is the direct-diagonalization result ||35|. The 
discrete transitions found in the direct diagonalization have been smeared with a gaussian 
of width 0.25 MeV in order to facilitate comparisons. While the SMMC total strength 
{i.e., the area under the curve) B{GT^) = 0.72 ± 0.11 [|T^ compares very well with the 
direct-diagonalization value of 0.79 |2^, the SMMC can recover only gross features of this 
distribution. In particular, the peak is somewhat too narrow, mainly due to the information 
lost by the Laplace transform. This attribution was checked by calculating the response 
function R{t) for the direct diagonalization distribution. (The peaks were smeared by 
Gaussians of 0.25 MeV width to account for the SMMC finite discretization.) This response 
function is shown in the lower left panel of Fig. 1, and agrees well with the SMMC result. 

The lower right-most panel in Fig. 1 shows the energy dependence of the cumulative 
strength, /g^ S{E')dE', where E* is the excitation energy in the daughter. One can see that 
the SMMC recovers the centroid and the width of the distribution reasonably. 

A brief remark about the possible sources of error is in order. Since our MaxEnt proce- 
dure provides a most probable extraction of the strength function, the strength distributions 
do not have error bars associated with them. However, from the SMMC error bars for R{t), 
we estimate the error in the position of the centroid to be about 0.5 MeV. In addition, 
we note that the response functions are measured in the parent nucleus, and to obtain the 
energy in the daughter we use the experimental mass excesses and a parametrization of the 



9 



Coulomb energy as defined in p3|. [In the test case (^^Ti), we exactly calculate this mass 



difference.] This parametrization provides a good overall description of the masses of the 



nuclei in this region |jT9[. We find an average deviation between 0.1 MeV (for A = A8 nuclei) 
and 0.5 MeV (for A = 54 nuclei) of our calculated binding energies from experimental values, 
suggesting that our procedure is quite justified. 

The upper panels of Fig. 1 show our results for the GT_ operator in ^'^Ti. The total 
strength, B{GT^), can be readily obtained from the renormalized Ikeda sum rule, B{GT^) — 
B{GT+) = 3(iV-Z)/(1.26)2 which is obeyed by both the SMMC and direct-diagonalization 
calculations. The GT_ operator takes the N > Z parent nucleus (with T = + 1) to 
T = (dotted), T = + I (dashed), and T = + 2 (not shown) states in the ^^V 
daughter. The T = states are the lowest in energy and contain most (85% in this case) 
of the strength. Assuming in the default model that the centroid of the T = + 1 states is 
located 5 MeV higher than the centroid of the T = Tz states, we obtain a good reproduction 
of both components of the strength distribution. This general assumption is experimentally 
valid in the even-even nuclei in this region. We also see at low energy a hint of the discrete 
low-energy states in the reconstruction. 

B. Comparison with experiment 

Experimental GT distributions are obtained from intermediate-energy charge exchange 
{n,p) [or {p,n)] cross sections at forward angles, which are proportional to the GT 
strength |3^. These experimental distributions typically extend only to 8 MeV in the 
daughter nucleus to exclude contributions from other multipolar it ies. 

We first compare our ^^Ti result for the GT+ distribution against experiment, as shown 
in Fig. 2. To simulate the finite experimental resolution and presentation of the data, 
the SMMC results have been smeared with Gaussians of standard deviation of 1.77 MeV, 



following Ref. ||3^. Our results are represented by the dotted line in Fig. 2, while the 



diagonalization results are shown as a solid histogram. The smeared diagonalization result 
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is shown by the dashed hne in the figure. The experimental B{GT^) distribution, shown 
as sohd dots, sums 1.42 ± 0.2 ^ compared to our renormahzed value of 0.71 ± 0.11. We 
find that the calculated Ohuj GT^ strength extends only over the region E* < 8 MeV (in 
agreement with the experimental value for this range of energy B{GT^) = 0.77 ± 0.1). This 
suggests that the GT+ strength observed for E* > 8 MeV corresponds to correlations outside 
our Ohuj model space. A similar conclusion has been reached in Ref. ||23|. We note that the 



inadequacy of a Ohu model space to describe the GT+ distribution at E* > 8 MeV might 
have some relevance to the f3f3 decay of ^'^Ca where considerable 2i//3/3 strength could 



be obtained from the overlap of this distribution with that of Ca in the (p, n) direction 
for these energies. However, the measured 21/(3(3 decay rate of ^^Ca p9i agrees well with the 



calculation based on the Ohu shell model, which includes the 1/1.26 normalization of the 
GT transition operator. 

We now turn to a comparison of SMMC results with experiment for nuclei in the mid-/p- 
shell where complete direct-diagonalization calculations are not possible. We first consider 
the (n, p) reaction and in Fig. 3 we show our results for all even-even nuclei with A = 48 — 64 
for which data are available The SMMC results have been smeared with Gaussians 

of standard deviation of 1.77 MeV to account for the finite experimental resolution, following 
Ref. [^. Experimentally, the GT+ strength is significantly fragmented over many states; 
the centroids and the widths of these distributions are reproduced very well in the SMMC 
approach. Our results for the total strengths are given in Table I. They agree with the 
data very well except for ^^Ni where our calculation underestimates the total experimental 
strength suggesting the need to augment the model space with the (79/2 and the (77/2 
orbitals. This shortcoming of the present model space is also visible in the GT+ distribution, 
which places the GT_|_ peak approximately 1.5 MeV above the experimental peak and misses 
the second peak at E^ 6 MeV, which is possibly due to ((79/2 — 97/2) transitions. 

SMMC results for odd-A nuclei in the {n,p) direction are shown in Fig. 4, where again 
the centroids and widths of the distributions are in good agreement with the data [|7|, ^ , P 



Calculations for odd-A nuclei are performed at a finite temperature of 0.8 MeV. (The tem- 
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perature dependence of these distributions will be discussed later in Section [IV Q .) The 
response functions for the three nuclei in Fig. 4 are sampled from the partition functions of 
their neighbors, i.e., ^^V from ^^Cr, ^^Mn from ^^Fe, and ^^Co from ^°Ni. The peaks of the 
observed GT+ distributions in odd-A nuclei in Fig. 4 are consistently at higher excitation 
energies in the daughter compared to the even-even cases in Fig. 3, a feature reproduced 
by the SMMC calculations. These higher excitation energies cause some Ohu strength to 
lie above the typical 8 MeV cut-off in odd-A nuclei. The data for ^^V and ^^Co have been 



analyzed for additional strength above 8 MeV p2| , p3i (see Table 1), while, to our knowledge, 
^^Mn has not been reanalyzed for potential GT strength at E* > 8 MeV. For even-even nu- 
clei the Ohuj GT_|_ strength appears to be exhausted at energies below 8 MeV, in agreement 
with the SMMC results shown in Fig. 3. Our results for ^^V and ^^Mn show some strength 
above 8 MeV, but this is not the case for ^^Co. 

In Fig. 5 we compare the GT_ distributions for a few nuclei where experimental data are 
available P,^; the experimental data for ^^Fe have been obtained from Ref. ||3^. From the 
cumulative strengths in the right panels of Fig. 5, we can conclude that the SMMC approach 
reproduces the experimental distribution moderately well for the cases of ^^Fe and ^^Fe. For 
the Ni isotopes, only partial information is available about these distributions. For ^^Ni the 
peaks in the experimental data |^ shown are to be associated with a finite width 1.3 MeV, 
0.7 MeV, and 0.5 MeV for the peaks at 9.2 MeV, 11.2 MeV, and 13.0 MeV, respectively. 
The strength in the giant resonance region between 6.4 MeV and 13.0 MeV is quoted as 5.5, 
while we obtain 6.1, which is consistent with the uncertainty in the excitation energy. For 
^°Ni the experimental value of the total GT_|_ strength [^] is 7.2 ± 1.8 whereas we obtain 
10.87 ± 0.23. As our calculation obeys the renormalized Ikeda sum-rule and reproduces the 
measured GT+ strength, the lower experimental value indicates some strength outside the 



experimental window of E* > 14 MeV. We also note that while Ref. |^ quotes an integrated 
strength of 6.22 between 4.0 and 14.0, MeV we obtain a value of 4.65. 
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C. Temperature dependence of GT strengths 



We now turn to the temperature evolution of GT+ strength functions. Representative 
strength distributions for two nuclei, ^^Co and ^''Ni, at several temperatures are shown in 
Fig. 6. Both figures are plotted as a function of E, the energy transfer to the parent nucleus. 
We note that the restriction of the model space to only /p-shell renders our calculation 
quantitatively unreliable for even-even nuclei at T ^ 1.4 MeV ||20|, while for the odd- A 
cases this temperature is likely even lower. 

With increasing temperature, three distinct effects occur that influence the GT strength 
distributions. 

• The number of states contributing to the thermal ensemble increases. Due to the 
pairing gap in even-even nuclei, this occurs at a higher temperatures in even-even 
nuclei than in odd-A nuclei. 

• GT transitions which are Pauli blocked at low temperatures due to closed neutron 
subshells (e.g. the /7/2 orbital) can become thermally unblocked as neutrons are moved 
to excited orbitals with increasing temperature. Similarly, protons which are thermally 
excited to higher orbitals can undergo allowed GT transitions. 

• The ground state in even-even nuclei is dominated by like-nucleon pairing. As indi- 
cated by SMMC calculations, these pairs break at around T = 1 MeV. Thus at low 
temperatures, a GT+ transition involves breaking a proton pair associated with an 
extra energy of 1-2 MeV. This "penalty energy" is removed at higher temperatures in 
states of higher excitation energy, in which the pair correlations are diminished. 

As we will discuss in the following, these three effects allow for an understanding of the 
temperature dependence of the GT_|_ strength distributions. 

In the case of ^^Co, with increasing temperature, the entire distribution shifts to lower 
excitation energies. The total strength decreases and the width of the distribution increases 
marginally with increasing temperatures. (We have checked that in the high-T limit, B{GT^) 
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rises to the single-particle value as expected.) Due to the lack of pairing of the odd particle 
in an odd-A nucleus, states of various spins are more quickly populated than in the even- 
even systems. These states then make transitions to daughter states by the GT operator. 
Thus, a plethora of states is easily accessible at moderate temperatures, and the required 
excitation energy in the daughter is lower. 

For ^''Ni, the peak in the strength distribution remains roughly constant with increasing 
temperature, while the width increases with the appearance of low-lying strength due to 
transitions from the thermally occupied to the empty excited orbitals. Note also that the 
centroid of the distribution remains constant at the low temperatures and then shifts to 
lower excitation at higher temperatures. The near constancy of the peak position in ^°Ni 
at low temperatures supports the shifting assumption [attributed to D. Brink in Ref. 



which states that the centroid corresponding to each parent excited state is shifted upward 
in the daughter nucleus by the energy of the parent state This hypothesis assumes 

that the internal configuration of the low-lying states is roughly the same. With increasing 
temperature, however, states with other internal configurations gain statistical weight, and 
in particular, the pair correlations in these excited states decrease. SMMC calculations 
indicate that pairs break around T = 1 MeV in even-even nuclei, allowing for a dramatic 
increase in thermally populated states in the parent at and above this temperature. For 
these excited states, no coherence energy has to be paid as penalty to break a proton pair in 
the GT transition, and the peak in the GT distribution moves to smaller energies. We also 
note that at temperatures T < 1.3 MeV the thermal ensemble already includes the lowest 
excited T + 1 states allowing for transitions at E = 0. In contrast, these transitions are not 
observed in ^^Co at the temperatures considered here, since the T + 1 states in this nuclei 
are at higher excited energies due to the larger neutron excess. We also observe a gradual 
decrease of the peak position with temperature in accordance with the fact that no pairing 
gap has to be overcome in odd-A nuclei. 
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V. SUMMARY AND CONCLUSIONS 



As mentioned in the Introduction, electron capture on iron region nuclei plays an impor- 
tant role at the onset of core collapse in a massive star. Under these conditions, nuclei have a 
finite temperature of 0.2—0.6 MeV. It is well known that for nuclei with an opened /p-shell 
neutron configuration, GT^ transitions dominate the electron capture rate, and a strong 
phase-space dependence makes the rate sensitive to the full GT+ distribution, rather than 
only to the total strength. Unfortunately, the GT+ strength is not experimentally accessible 
for those nuclei of importance in the presupernova collapse. Thus, collapse studies have to 
rely on theoretical estimates which, until recently, could not be performed with great con- 
fidence. This has now changed. As SMMC calculations reproduce the measured data from 
first principles without nucleus-specific data fitting (which has been necessary in previous 
studies), they are reliable enough to predict the GT+ distributions for those astrophysically 
important nuclei not experimentally accessible. SMMC calculations for these nuclei are in 
progress. 

In this paper, we have calculated response functions for the Gamow- Teller operators for 
several nuclei in the /p-shell. We use the KB3 interaction, which is well suited for Ohu 
calculations. Using an implementation of the MaxEnt technique, we have then obtained the 
corresponding strength distributions. 

The extracted Gamow- Teller distributions compare very well with both direct- 
diagonalization calculations and the experimentally obtained distributions. We note that we 
invoke the standard renormalization factor of 1/1.26 for the transition operator, in keeping 
with the observation in sd— and /p-shell nuclei that complete Ohu calculations require this 
overall renormalization for agreement with experiment. 

We have also studied the effect of finite temperature on Gamow- Teller distributions 
and have demonstrated for sample nuclei that our calculations at T = 0.8 MeV should be 
adequate to describe the distributions required to calculate electron capture rates for the 
pre-supernova problem . Studies of the Gamow- Teller strengths and electron capture rates 
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for nuclei relevant to the presupernova collapse will be described elsewhere. 
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TABLES 

TABLE L Renormalized B{GT^) strengths as calculated in the SMMC approach compared to 
experimental strengths P,^, [40| , ^ . The superscripts on the experimental results indicate the upper 
limit of energies used to obtain the total strength, a: up to 14 MeV. b: up to 12.5 MeV. c: up to 
10 MeV. d: up to 8.5 MeV e: up to 8 MeV 



nucleus 


fi(GT+)(SMMC) 


fi(Gr+)(expt) 




0.71 ±0.11 


1.31 ±0.2" 


Sly 


1.40 ±0.14 


1.48 ±0.03^ 


54pg 


3.84 ±0.28 


3..1 ±0.6" 


55Mn 


1.84 ±0.36 


1.7 ±0.2^ 




2.51 ±0.17 


2.9 ±0.3"^ 


58Ni 


4.23 ±0.31 


3.8 ±0.4^ 




2.60 ±0.31 


2.39 ±0.07^ 


60Ni 


3.26 ±0.25 


3.11 ±0.08" 


62Ni 


2.16 ±0.25 


2.53 ±0.07" 


64Ni 


1.09 ±0.18 


1.72 ±0.09" 
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FIGURES 

FIG. 1. Left-most panels show the GT_ (upper) and GT_(_ (lower) response functions cal- 
culated through the SMMC. The middle panels show the corresponding strength function and 



direct-diagonalization results [^,35] in the corresponding daughter. For the GT_ we show both 
the T = Tz and T = Tj + 1 channels, while the dash-dot line in the GT+ distribution comes 
from folding the SMMC results with a gaussian corresponding to the experimental sensitivity. The 
right-most panels show the cumulative strengths as a function of the daughter excitation energy. 
For the GT_ we show the cumulative T = + 1 strength starting from the total in the T = Tz 
channel. 

FIG. 2. Calculated strength function (smeared by the experimental resolution) for the GT+ 
operator for ^^Ti compared to the experimental data H]. Also shown is the shell model spectrum 
obtained by diagonalization, and smeared by 0.25 MeV (histogram) and by 1.77 MeV (dashed line) 
to account for the experimental resolution. 



FIG. 3. Comparison of calculated GT+ strength distribution against experiment |Q,^,40,41| 
for even-even nuclei as function of excitation energy in the corresponding daughter nuclei. 

FIG. 4. Calculated GT+ distributions for odd-^ nuclei. Also shown are the experimental 
distributions [0g|,|l|. The energies are in the corresponding daughter. 



FIG. 5. Left panels: Calculated GT_ distributions for several nuclei in the mid-/p shell against 
distributions obtained from {p,n) reactions []^J4^. Right panels: Cumulative strength distributions 
versus daughter excitation energy for SMMC calculations and experiment. B{GT-) from SMMC 
(solid circles) and from experiment (open circles) are shown staggered for clarity. 

FIG. 6. Temperature evolution of GT-|- strength distribution for sample nuclei (left: ^^Co; 
right: ^^Ni) versus parent excitation energy. 
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